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1 .  INTRODUCTION 


Surface  roughness  can  play  an  important  role  in  increasing  friction 
and  heat  transfer  under  turbulent  boundary  layer  conditions.  Roughness  ef¬ 
fects  have  been  studied  extensively  over  the  past  fifty  years,  in  connection 
with  applications  such  as  head  losses  in  pipes  and  other  hydraulic  equipment, 
ship  hull  drag,  airplane  drag,  and  high  speed  missile  drag  and  heating.  This 
study  emphasizes  the  effects  of  roughness  character  euid  the  behavior  of  rough- 
wall  boundary  layers  under  supersonic  and  hypersonic  conditions.  Hypersonic 
boundary  layers  tend  to  be  quite  thin,  which  naturally  increases  the  likeli¬ 
hood  that  roughness  will  be  important.  There  have  been  few  experimental  or 
theoretical  investigations  of  rough  wall  turbulent  boundary  layers  at  high 
Mach  numbers.  Furthermore,  many  high  speed  flight  vehicles,  such  as  re-entry 
vehicles,  are  fabricated  from  composite  materials.  Since  these  typically  in¬ 
volve  woven  fibers  filled  with  a  resin,  they  present  a  different  roughness 
character  than  would  be  found  on,  say,  a  metallic  surface. 

One  cannot  adequately  discuss  the  behavior  of  flows  over  rough  sur¬ 
faces  without  recalling  the  classic  experiments  of  Nikuradse,^  in  which  water 
was  flowed  through  pipes  roughened  by  sand.  In  the  "fully  rough"  regime,  the 
measured  friction  factor  X  was  found  to  depend  only  on  the  ratio  of  roughness 
height  kg  and  pipe  radius 

X  =  [1.74+2  logio  R/ks]“2  .  (1) 

(A  list  of  symbol  definitions  is  given  on  page  61).  This  result  holds  for  kg'*" 
=  U-rkg/v  greater  than  about  70.  Smooth  wall  behavior  prevails  for  kg'*’<  5,  and 
Nikuradse  presented  a  correlation  for  the  intermediate  transition  regime. 
Analogous  experiments  were  performed  by  Dipprey  and  Sabersky^  to  extend 
Nikuradse 's  results  to  heat  transfer  at  various  Prandtl  numbers.  Although 
Nikuradse 's  sand  grains  were  carefully  sifted  to  obtain  a  relatively  uniform 
size  (diameter  =  kg),  we  have  no  detailed  information  on  the  statistics  of  the 
roughened  surface  that  resulted  from  adhesively  bonding  the  sand  grains  to  the 
smooth  surface. 

Some  very  detailed  measurements  were  obtained  by  Moffat  and  co- 
workers^"®  on  the  low  speed  flow  of  air  over  a  flat  plate  covered  by  closely 


packed  spheres.  Their  data  Include  skin  friction,  heat  transfer,  and  profiles 
of  mean  and  fluctuating  quantities.  These  results  are  quite  useful  for  vali¬ 
dating  theoretical  models.  However,  only  one  roughness  was  investigated. 

Another  important  class  of  experiments  involves  "two-dimensional" 
roughness,  such  as  machined  grooves  or  square  rods,  normal  to  the  flow  direc¬ 
tion.  Betterman®  varied  the  relative  rod  spacing  by  a  factor  of  about  three, 
and  Antonia  and  Luxton^  made  detailed  surveys  of  the  turbulence  parameters  in 
the  boundary  layer  over  this  type  of  roughness.  While  several  authors  have 
correlated  2-D  roughness  data  along  with  3-D  or  distributed  roughness  data, 

2-D  roughness  will  not  be  considered  here.  We  might  expect  substantial  dif¬ 
ferences  in  the  nature  of  the  flows.  With  2-D  roughness,  the  flow  is  likely 
to  be  dominated  by  cavities  in  the  grooves  between  elements,  whereas  separa¬ 
tion  should  be.  less  important  with  distributed  roughness.  The  model  presented 
below  is  aimed  entirely  at  distributed  roughness,  appropriate  to  the  vast 
majority  of  practical  applications. 

The  available  data  base  on  roughness  effects  in  compressible  flows  is 
considerably  smaller  than  that  for  low  speeds.  The  most  extensive  set  of  ex¬ 
periments  were  those  of  the  Passive  Nosetip  Technology  (PANT)  program,®  in 
which  roughened  hemispherical  models  were  placed  in  NSWC  Tunnel  8  at  a  free- 
stream  Mach  number  Mo,  =  5.  Heat  transfer  was  measured  by  calorimeter  methods. 

The  roughness,  which  varied  over  two  orders  of  magnitude  in  height,  was  creat¬ 
ed  by  grit  blasting  or  by  bonding  grit  particles?  a  considerable  quantity  of 
data  were  obtained  on  both  roughness  augmentation  of  turbulent  heating  as  well 
as  roughness-induced  transition.  It  should  be  emphasized  that  these  tests 
were  performed  on  blunt  nose  regions,  where  the  boundary  layer  edge  Mach  num¬ 
bers  are  subsonic  or  modestly  supersonic.  While  there  is  a  substantial  varia¬ 
tion  of  density  and  temperature  across  the  boundary  layer  in  these  conditions, 
they  certainly  are  not  representative  of  high  Mach  number  boundary  layers. 

Holden®  has  run  tests  in  the  hypersonic  shock  tunnel  at  Calspan 
(M«b  =  11-13)  on  45°  cones  to  which  grit  was  bonded;  the  boundary  layer  edge 
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Mach  number  is  about  1.8.  Higher  edge  Mach  numbers,  about  4.8,  were  obtained 
in  NSWC  Tunnel  2  by  Keel,^®  on  5°  cones  with  sand  grains  attached  by  epoxy. 


In  the  S2une  facility^  Volslnet^^  measured  the  combined  effects  of  roughness 
and  mass  addition.  He  created  the  roughness  by  covering  a  porous  section  of 
the  tunnel  wall  with  a  screen.  This  provides  a  questionable  simulation  of 
distributed  roughness,  but  is  undoubtedly  necessitated  by  the  extreme  diffi¬ 
culty  in  fabricating  a  rough  porous  surface  with  controlled  permeability. 

Truly  hypersonic  tests  have  been  performed  by  Hill (nswc  Hyperveloc¬ 
ity  Tunnel)  and  Holden^^  (Calspan  hypersonic  shock  tunnel),  using  grit  bonded 
to  the  surface  of  slender  cones.  Hill's  experiments,  on  a  7®  cone  at  Mg  * 

8.1,  used  three  different  roughness  heights.  Holden  employed  a  single  rough¬ 
ness  on  a  6®  cone  but  introduced  varying  angles  of  attack  to  thin  the  boundary 
layer  on  the  windward  side,  effectively  increasing  the  roughness  effect  (Mg  « 
4.4  -  9.4).  Despite  the  general  similarity  of  the  conditions  for  these  two 
experii.ients,  there  are  differences  in  thu  results  that  deserve  further  discus¬ 
sion  below. 

Investigations  on  the  effects  of  roughness  character,  where  the  shape 
and  spacing  of  roughness  elements  are  varied,  are  rather  limited.  "Hie  classic 
experiment  was  performed  by  Schlichting, on  one  wall  of  a  water  channel. 
Various  arrangements  of  roughness  elements  were  used.  Including  spheres, 
spherical  segments,  cones,  and  short  angles,  at  several  relative  spaclngs. 
Results  were  presented  in  terms  of  the  wall  shear  and  equivalent  sand  grain 
roughnesses.  Several  others  have  investigated  the  roughness  character  effect 
over  more  limited  ranges,  mostly  by  changing  the  spacing  for  a  given  shape. 
Chen  and  Roberson studied  three  relative  spacings  of  hemispheres,  for  air 
flow  through  a  pipe.  Raupach,  Thom  and  Edwards^®  varied  the  relative  spacing 
over  a  factor  of  four  for  cylindrical  elements  (k/D  =  1)  in  a  wind  tunnel.  In 
water  channel  flume  experiments,  Mlrajgaoker  and  Charlu^^  used  stones  at  six 
spacings  and  Sayre  and  Albertson^^  used  baffles  (similar  to  Schlichting' s 
short  angles,  but  with  width  =  4  x  height),  again  with  six  spacings.  O'Lough- 
lin  and  Annambhotla^^  used  cubes  at  three  spacings,  in  a  wind  tunnel.  Mul- 
hearn  and  Finnigan^®  examined  gravel  at  one  spacing  in  a  wind  tunnel.  To  sim¬ 
ulate  plant  crop  canopy  flow,  Thom^^  and  Seginer  et  al.^^  studied  tall,  thin 
rods  (k/D  -  100),  which  provide  a  major  variation  in  aspect  ratio  of  the 
roughness  elements. 


One  series  of  roughness  character  tests  under  supersonic  conditions 
has  recently  been  conducted  by  Acurex  Corp.^^  in  AEDC  Tunnel  F.  The  models 
were  45°  cones,  with  an  edge  Mach  number  of  1.7.  Seven  surfaces  were  used: 
essentially  smooth,  grit  blasted,  bonded  grit,  and  four  chemically-etched 
roughness  patterns  (two  heights,  two  spacings  each).  Ibe  accuracy  of  the 
measured  heating  rates  may  be  limited  by  the  fact  that  Tunnel  F  was  an  arc 
heated,  hot-shot  type,  with  pressure  decreasing  continuously  during  the  test 
time,  and  by  the  fact  that  the  roughness  characteristics  varied  over  the  model 
surface. 

A  great  number  of  methods  of  varying  degrees  of  empiricism  have  been 
developed  to  predict  rough  wall  drag  and  heating.  The  soundest  methods  start 
with  Nikuradse's  observation  that  the  logarithmic  portion  of  the  mean  velocity 
is  shifted  downward  by  roughness.  This  downward  shift,  AU.j/U^,  is  also  di¬ 
rectly  related  to  the  roughness-induced  increase  in  friction  by  the  following 
relation  (at  least  for  low  speed  flows) 


The  velocity  shift  depends  on  wall  conditions,  and  thus  is  a  function  only  of 
k"*".  The  smooth  wall  friction  coefficient  generally  depends  on  the  Reynolds 
number  based  on  some  measure  of  the  thickness  of  the  viscous  zone  (e.g.  Re0), 
as  well  as  on  pressure  gradient  and  other  geometrical  factors  (boundary  layer 
vs.  pipe,  etc.).  In  the  fully  rough  regime,  however,  the  dependence  of  C 

^sm 

on  Re0  combines  with  dependence  of  AUi/U^  on  k"*"  to  yield  a  dependence  only  on 
the  ratio  k/9  (or  k/R  in  a  pipe),  independent  of  Reynolds  number.  Note  that 
Eq.  (2)  is  transcendental  in  C^,  since  k"*"  involves  Also,  the 

above  equation  cannot  be  applied  directly  to  compressible  flows. 

For  sand  grain  roughness,  the  velocity  shift  in  the  fully  rough  regime 
is 


5.6  log  kg 


3 


(3) 


U 

T 


Dvorak24»25  first  extended  this  relation  to  describe  variations  in  roughness 
character 


AU 

^ —  »  5.6  log  k'*’  +  f{X)  ,  (4) 


so  that  the  equivalent  sand  grain  roughness  (kg)  is  related  to  the  actual 
roughness  height  (k)  by 


5.6  log  kgA  =  f(  X)  +  3  .  (5) 

Here.  f(X)  is  some  measure  of  the  roughness  density.  Dvorak  correlated  the 
available  data  with  X“^  being  the  fraction  of  the  surface  covered  by  roughness 
(the  base  area  of  roughness  elements  per  unit  underlying  smooth  surface  area). 
Simpson^®  obtained  improved  results  in  terms  of  X)^,  the  inverse  of  the  pro¬ 
jected  frontal  area  of  the  elements  per  unit  surface  area,  and  others^'^'^® 
have  used  similar  variations  of  the  Dvorak  approach. 

Figure  1  shows  Simpson's  roughness  density  correlations.  The  two 
straight  line  segments  are  the  same  as  Dvorak's  (in  terms  of  X  rather  than 
Xjj).  Note  that  the  roughness  effect  has  a  maximum  at  =  4.70,  which  corres¬ 
ponds  to  spherical  elements  with  an  average  separation  of  1.9  diameters.  Note 
also  that  nearly  all  of  the  data  to  the  left  of  this  maximum,  for  more  closely 
packed  elements,  are  for  two-dimensional  roughness.  Simpson  notes  that  the 
flow  in  this  regime  may  be  dominated  by  cavity  flow  between  the  elements,  and 
speculates  that  this  decrease  in  roughness  effect  with  decreasing  spacing  may 
not  occur  for  3-D  roughness.  We  shall  return  to  this  matter  in  some  depth  be¬ 
low. 

For  compressible  flows,  Dvorak^®  applied  the  same  compressibility  fac¬ 
tor  used  for  smooth  walls  (these  are  surveyed  in  Ref.  29).  Equation  (2)  is 
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.k. 


used  to  compute  an  incompressible  friction  coefficient/  with  ^fgg,  based  on 
F0*Re0  and  AUi/U-j  based  on  the  actual  wall  conditions;  the  compressible  fric¬ 
tion  coefficient  is  then  obtained  by  dividing  the  incompressible  value  by  Fq. 
However,  this  procedure  is  not  necessarily  logically  consistent  with  the  in¬ 
compressible  behavior.  If  AUi/U^  really  depends  on  local  wall  conditions, 
then  Eg.  (2)  would  apply  directly  to  compressible  situations,  simply  by  in¬ 
serting  the  appropriate  compressible  value  of  ^fsm*  Dvorak  found  this  al¬ 

ternative  to  yield  poor  agreement  with  data. 

A  number  of  investigators  have  developed  more  empirical  methods  for 
reentry  applications.  One  of  these,  which  is  well-known  and  has  been  derived 
from  both  low  speed  and  high  speed  measurements,  is  the  version  developed  by 
Acurex  Corp.  emd  contained  in  the  ABRES  Shape  Change  Code.^^  The  roughness 
augmentation  of  skin  friction  is  described  by  an  influence  coefficient,  by 
which  the  smooth  wall  value  should  be  multiplied,  independent  of  compressibil¬ 
ity; 

If  =  1  +  0.5  fi(k/e)gi(x)  (6) 

fl(k/e)  *  1  +  0.09  k/e  +  0.53(  1-e"’'/®) 

91 'X)  “  X  +  1.5(1-e"X)  for  x  ^  0 

»  0  for  X  ^  0 

X  =  log (k'*’/ IS. 5)  (Based  on  smooth  wall  Cf). 

While  this  relation  provides  a  smooth  transition  between  smooth  and  rough  wall 
behavior,  it  does  not  obviously  reduce  to  the  dependence  observed  by  Nikuradse 
in  the  fully  rough  regime,  in  that  the  drag  is  predicted  to  depend  on  both  k"*" 
and  k/e.  Dahm^^  has  recently  developed  a  new  correlation  for  roughness  ef¬ 
fects,  involving  a  "two-layer"  model. 

Another  issue  of  considerable  uncertainty  is  the  relation  between  heat 
transfer  and  skin  friction.  For  smooth  walls,  the  mechanisms  of  turbulent 
heat  emd  momentum  transfer  are  quite  similar  and  the  heat  transfer  and  skin 
friction  coefficients  are  closely  related.  But  with  rough  walls,  form  drag 
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on  the  elements  has  no  analogous  thermal  mechanism.  This  provides  an  intui¬ 
tive  explanation  for  the  common  observation  that  roughness  causes  smaller  In¬ 
creases  in  heating  than  in  friction.  Owen  and  Thompson^^  developed  a  correla¬ 
tion  for  roughness-dominated  heat  transfer,  following  lines  suggested  by  as¬ 
suming  cavity  flow  between  the  elements.  However,  the  result  involves  prima¬ 
rily  k"*"  and  it  is  not  clear  whether  this  behavior  is  consistent  with  the  well 
known  dependence  of  skin  friction  on  k/6  rather  than  k"*".  The  heat  transfer 
correlation  used  in  the  ABRES  Shape  Change  Code^®  is  simple  and  pragmatic  -  the 
numerical  factor  0.5  in  Eq.  (6)  is  replaced  by  0.3. 

The  analyses  to  be  presented  below  are  based  on  a  fairly  basic  model 
for  rough  wall  boundary  layers,  wherein  a  Reynolds  stress  turbulence  model  is 
combined  with  a  form  drag  'description  for  the  effect  of  roughness  elements  on 
thii  flow.  The  computer  model  is  used  ir.  tv/o  ways;  1)  it  is  compared  against 
relevant  data  to  establish  the  validity  and  accuracy  of  the  theory  and  to  of¬ 
fer  explanations  for  observed  trends;  and  2)  it  is  used  to  develop  correla¬ 
tions  for  engineering  applications  by  suggesting  scaling  laws  and  by  providing 
nimierical  data  to  be  correlated.  In  previous  papers, we  showed  con5>ari- 
sons  with  a  good  portion  of  the  available  measurements  and  presented  some  pre¬ 
liminary  scaling  laws.  Here  we  will  analyze  some  of  the  most  important  data 
on  roughness  character  and  conpressibility,  and  shall  present  correlations 
that  should  have  widespread  utility. 
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2.  ROUGH  WALL  BOUNDARY  LAYER  MODEL 


The  basic  model  for  rough  wall  turbulent  boundary  layers  is  the  same 
as  that  used  previously, and  we  shall  only  outline  the  most  relevant  fea¬ 
tures  here.  The  turbulence  model  is  a  Reynolds  stress  or  second-order  closure 
method,  which  confutes  both  mean  and  fluctuating  velocities  and  temperatures. 
The  dependent  velocity  variables  are  the  mean  velocity  vector  Uj^,  the  Reynolds 
'tress  tensor  Uj^'uj',  and  the  isotropic  dissipation  rate  4.  The  analogous 
thermal  variables  (temperature  or,  more  precisely,  enthalpy  h)  are  the  mean 

enthalpy  h,  the  mean  square  fluctuating  enthalpy  h'^,  and  the  Reynolds  heat 
flux  vector  u^'h*.  Under  the  boundary  layer  approximation,  this  set  of  vari¬ 
ables  reduces  to  U,  V,  u^,  v'^,  w*^,  u'v',  4,  h,  h'^,  u'h',  and  v'h*.  The 
closure  approximations  that  are  used  to  derive  the  required  equations  are 
somewhat  standard  at  this  time.  The  formulation  has  been  successfully  applied 
to  a  variety  of  smooth  wall  boundary  layer  and  free  shear  flows.  Predicted 
smooth  wall  skin  friction  coefficients  are  generally  within  10-15%  of  accepted 
values,  except  for  a  few  cases  with  very  large  density  differences,  to  be 
noted  below. 

Our  model  makes  the  basic  assumption  that  the  force  on  roughness  ele¬ 
ments  can  be  viewed  as  form  drag.  This  implicitly  requires  that  the  flow  ap¬ 
proaching  an  individual  element  be  attached.  As  already  noted,  cavity  flow  is 
likely  to  prevail  with  2-D  roughness;  for  this  reason  the  present  model  should 
be  more  appropriate  for  distributed  roughness. 

The  rough  surface  is  idealized  as  being  made  up  of  identical  elements 
(although  the  extension  to  a  size  distribution  would  be  straightforward).  The 
bottom  of  the  elements,  or  the  underlying  smooth  wall,  is  at  y=0.  The  element 
height  is  k,  and  I  is  the  average  element  spacing  (f”2  gives  the  number  of 
elements  per  unit  area) .  We  restrict  our  treatment  to  elements  with  circular 
cross  sections  at  all  heights,  with  D(y)  denoting  the  diameter  at  height  y  (y 
<  k) .  As  discussed  in  more  detail  in  Ref.  33,  form  drag  on  the  elements  is 
described  by  an  appropriate  negative  (sink)  term  in  the  mean  momentum  equa¬ 
tion: 


Ru  =  -  pu2Cd  D(y)/«2  . 


(7) 
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A  drag  coefficient  value  of  Cq  =  0.6  is  roughly  appropriate  for  elements  such 
as  cones  or  hemispheres.  In  addition,  there  should  be  source  terms  for  turbu¬ 
lent  kinetic  energy  and  dissipation,  describing  the  tendency  of  roughness  to 
increase  velocity  fluctuations.  These  terms,  which  are  discussed  in  Ref.  33, 
are  not  very  important  compared  to  the  indirect  effect  of  roughness  to  in¬ 
crease  the  turbulent  energy  by  increasing  the  mean  shear. 

Except  in  the  Stokes  flow  regime,  heat  transfer  to  an  element  should 
be  small.  Therefore,  the  only  roughness  term  appearing  in  the  thermal  equa¬ 
tions  is  a  source  term  for  the  mean  static  enthalpy.  "Ihis  term  is  constructed 
so  that,  in  combination  with  Eq.  (7),  form  drag  does  not  alter  total  enthal¬ 
py. 


Rh  =  +  1/2pu3Ci3D(y)A2  .  (8) 

It  is  also  necessary  to  account  for  the  blockage  effect  of  the  rough¬ 
ness  elements.  At  a  given  height,  the  fraction  of  the  flow  area  normal  to  the 
X  direction  is  1-D(y)/f.  Terms  that  act  in  the  streamwise  direction,  such  as 
the  convective  operator  pUd/dx,  are  multiplied  by  this  factor.  Terms  that  act 
on  planes  normal  to  the  y  direction,  or  that  act  on  a  unit  volume,  should  be 
modified  by  However,  the  loughness  terms  disrusred  above  aro  al¬ 

ready  based  on  the  total  volume,  rather  than  the  available  flow  volume,  and 
need  no  such  factor.  If  the  entire  equation  is  divided  by  B(y)  =  1-irD2/4Jt,2,  a 
relatively  simple  result  is  obtained.  For  excunple,  the  mean  momentum  equation 
becomes 

9U  3U  1  3  du 

f(y)pu  —  +  pv  —  •=  -  f  (y)  ^  +  -  i  Bjj  — 

3x  3y  3x  B  3y  3y 

-  —  (pu'v’)  -  -  pU^Cp  ^  B-’  (9) 

3y  2 
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where 


,  _  1-D/t 

1-irD^/4£^ 


(10) 


The  function  f(y)  contains  the  main  effect  of  bloc)cage  and  may  be  ab¬ 
sorbed  in  the  definition  of  the  stre£un  function  that  is  introduced  to  elimi¬ 
nate  the  normal  velocity 

g  =  f(y)pu  ,  |1  =  -  pv  .  (11) 

Note  that  if  the  elements  are  packed  so  tightly  that  they  are  touching  over 
some  range  of  y,  then  D  =  £  and  f(y)  =  0  over  that  range.  Our  formulation 
forces  the  velocity  to  remain  zero  up  to  the  height  where  D  <  £  and  the  flow 
is  unblocked.  Of  course/  common  sense  would  dictate  redefining  y  •  0  as  the 
lowest  point  where  the  flow  is  unblocked. 

A  major  advantage  of  this  model  is  that  solutions  are  obtained  for 
both  velocity  and  thermal  variables.  Heat  transfer  is  obtained  directly, 
without  invoking  a  Reynolds  analogy.  Finite  difference  solutions  are  obtained 
using  the  obvious  boundary  conditions.  Fluctuating  quantities  are  zero  at  the 
base  of  the  wall,  y  =  0.  At  the  outer  edge,  fluctuating  quantities  are  zero 
for  a  boundary  layer  or  obey  a  symmetry  condition  for  a  channel  flow.  For 
numerical  solutions,  the  equations  are  first  transformed  to  the  stream  func¬ 
tion  coordinate,  guaranteeing  mass  conservation  and  eliminating  the  normal 
velocity,  V. 

The  transverse  coordinate  is  normalized  by  the  edge  value  of  the 
stream  function  so  that  additional  mesh  points  need  not  be  carried  in  the  free 
stream  to  allow  for  boundary  layer  growth.  For  proper  resolution  of  the 
region  near  the  wall,  a  linear  mesh  in  the  logarithm  of  the  stream  function  is 
used.  The  finite-difference  equations  are  solved  with  a  block  tridiagonal 
Newton-Raphson  technique. 
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It  should  be  noted  that  the  use  of  such  a  model  is  not  unique  to  this 
study.  Lin  and  Bywater^S  have  used  the  present  form  drag  treatment  in  a  two- 
equation  TKE  model,  and  generally  obtained  results  superior  to  those  obtained 
from  lower  level  approaches.  A  very  similar  model  has  been  developed  by  Wil¬ 
son  and  Shaw^^  to  analyze  transport  processes  in  plant  canopies.  A  crop  of 
corn  plants,  or  a  stand  of  trees,  represents  an  attractive  geometry  for  the 
form  drag  model.  For  such  applications,  it  may  be  necessary  to  account  for 
the  effect  of  plant  motion  on  the  turbulent  energy  budget. 


3.  SPECIFICATION  OF  ROUGHNESS  CHARACTERISTICS 

An  important  aspect  of  predicting  or  analyzing  roughness  effects  is 
the  proper  specification  of  the  roughness  characteristics.  For  laboratory  ex¬ 
periments  in  which  identical  elements  of  a  simple  shape  (spheres,  cones,  etc.) 
are  attached  to  a  smooth  surface  in  a  regular  pattern,  the  required  element 
height,  shape  and  average  spacing  are  obvious.  Such  is  not  necessarily  the 
case  for  grains  of  sand  or  grit,  sifted  to  a  narrow  size  range  and  applied  to 
the  surface  with  an  adhesive.  Ihe  details  of  the  bonding  technique  and  the 
departure  from  sphericity  of  the  particles  can  affect  the  roughness  pareim- 
eters.  Real  surface  materials  introduce  even  greater  uncertainties.  Here  we 
present  the  method  for  deriving  the  roughness  specifications,  required  for  our 
model,  from  profiloraeter  surface  measurements. 

A  profilometer  measurement  of  the  surface  in  question  is  required. 

This  consists  of  an  irregular  trace  of  height  y,  above  some  reference,  as  a 
function  of  distance  along  the  line  over  which  the  stylus  was  traced.  It  will 
be  assumed  here  that  all  elements  are  identical  in  size  and  shape  -  one  could 
derive  a  more  sophisticated  analysis  for  situations  where  a  significant  varia¬ 
tion  in  element  sizes  is  expected.  It  is  also  assumed  that  location  of  rough¬ 
ness  elements  has  at  least  a  moderate  degree  of  randomness. 

Figure  2  sketches  a  typical  profilometer  curve.  From  this  curve,  it 
is  necessary  to  form  a  probability  of  exceedance  distribution,  P.E.(y),  which 
is  simply  the  fraction  of  the  trace  with  heights  greater  than  y.  One  must  al¬ 
so  identify  peaks  and  conpute  the  average  spacing  between  peaks,  L^.  Finally, 
one  must  define  the  height  y  =  0,  corresponding  to  the  effective  f’oor  of  the 
roughness. 

Now,  any  element  located  within  a  roughness  element  base  radius  on 
either  side  of  the  profilometer  line  should  be  detected  and  counted  as  a  peak. 
Since  is  the  number  of  elements  per  unit  area,  a  profilometer  trace  of 
length  should  detect  peaks.  The  number  detected  is  L^/1^,  by 

definition  of  the  average  peak  spacing.  Equating  these  two  quantities  gives 

=  D(0)Lp  .  (12) 
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Since  the  roughness  statistics  are  taken  to  be  uniform,  the  probabil¬ 
ity  of  exceedance  of  height  y  along  a  profilometer  trace  must  equal  the  frac¬ 
tion  of  the  surface  area  above  that  height.  The  latter  fraction  is  simply  the 
cross-sectional  area  of  an  element  at  the  height,  Tr/4  D^(y),  multiplied  by  the 
density  of  elements, 

P.E.(y)  =  J  D2(y)/£2  (13) 

Equations  (12)  and  (13)  are  sufficient  to  reconstruct  the  roughness 
characteristics.  From  the  P.  E.  value  at  y  =  0  (which  should  be  less  than 
Tr/4  to  allow  D(0)  <  t),  Eqs.  (12)  and  (13)  yield  D(0)  and  i.  Equation  (13) 
can  then  be  used  to  obtain  D(y)  for  greater  heights. 

Some  examples  of  the  roughness  specifications  that  result  from  this 
process  will  be  given  below,  for  grit-bonded  surfaces.  In  practice,  some 
judgement  is  required  to  identify  the  "floor"  of  the  surface  (y  =  0)  as  well 
as  to  identify  peaks  in  the  profilometer  curves.  With  sufficient  statistics 
(e.g.,  the  profilometer  traces  should  intersect  on  the  order  of  100  or  more 
elements),  the  practical  effects  of  such  uncertainties  is  small,  because  the 
portions  of  the  elements  near  the  base  do  not  contribute  greatly  to  the  total 
drag. 
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4.  COMPARISONS  WITH  ROUGHNESS  CHARACTER  DATA 

As  indicated  in  the  Introduction,  the  most  extensive  experiment  on 
roughness  character  is  the  Schlichting^**  study  of  various  roughness  patterns 
on  one  wall  of  a  water  channel.  In  Ref.  34  we  made  some  preliminary  analyses 
of  Schlichting' s  data,  using  boundary  layer  solutions  at  the  proper  value  of 
Re{.  Here  in  Figs.  3-6,  we  show  more  appropriate  solutions  for  the  proper 
channel  geometry. 

For  spheres.  Fig.  3,  the  rough  wall  boundary  layer  model  agrees  ac¬ 
ceptably  with  the  measured  friction  coefficients,  the  worst  error  being  about 
25%  at  =  3  or  l/D  ~  1.5.  At  the  closest  packing,  we  specified  hemispheres, 
since  the  bottom  half  of  spheres  would  be  completely  blocked.  This  is  done 
for  numerical  convenience,  to  avoid  singularities  associated  with  f(y)  =  0  in 
Eq.  9.  According  to  our  calculations,  blockage  of  the  bottom  half-sphere  is 
responsible  for  most  of  observed  reduction  in  drag  as  l/D  1 .  For  compari¬ 
son,  the  Dvorak/Simpson  correlations  (which  are  identical  for  spheres)  are 
also  plotted. 

For  spherical  segments,  which  are  somewhat  less  than  hemispheres,  pro¬ 
bably  to  simulate  rivets,  the  computer  code  gives  excellent  results  (Fig.  4). 
Note  that  the  code  shows  only  a  very  modest  maximum  in  Cf,  in  contrast  to  the 
Dvorak  and  Simpson  correlations.  Good  agreement  was  also  obtained  with  Chen 
and  Roberson 's^^  data  on  hemispheres,  although  they  investigated  only  very 
large  roughness  spacings. 

For  Schlichting' s  cones,  shown  in  Fig.  5,  our  model  is  slightly  low. 
Closely  packed  cones  were  not  tested,  but  again  the  present  model  shows  no 
pronounced  p>eak.  The  short  angles  of  Fig.  6  were  simulated  as  cylinders  of 
the  same  height  and  width  in  our  calculations.  The  model  is  noticeably  low 
here,  for  unknown  reasons.  However,  we  would  expect  less  accuracy  for  cases 
with  non-circular  roughness,  and  similar  errors  would  be  expected  for  the  baf¬ 
fles  of  Sayre  and  Albertson.^®  Generally  good  agreement  was  achieved  with  the 
data  of  Raupach,  Thom  and  Edwards,^®  on  cylinders  (height  i  diameter),  as 
shown  in  Fig.  7. 

Figure  8  compares  our  results  with  the  observations  of  Mirajgaoker  and 
Charlu. We  used  the  channel  flow  version  of  our  model  to  simulate  their 
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flume,  and  specified  their  stones  as  spheres.  It  is  obvious  that  the  theory 
significantly  underpredicts  the  measured  drag,  as  does  Sinpson's  correla¬ 
tion.  ^6  The  reasons  for  this  discrepancy  are  unclear  -  perhaps  free  surface 
effects  are  important.  However,  no  attenpt  has  been  made  to  investigate  that 
possibility  in  this  study. 
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5.  COMPARISONS  WITH  COMPRESSIBLE  ROUGHNESS  DATA 

The  effect  of  surface  roughness  in  conpressible  flow  conditions  re¬ 
mains  poorly  understood,  perhaps  largely  because  the  available  data  base  is 
rather  fragmentary.  In  Ref.  33  we  showed  extensive  conparisons  with  the  PANT 
data  on  hemispherical  nosetips  at  supersonic  velocities,  and  we  also  demon¬ 
strated^^  satisfactory  agreement  with  Keel's  data^®  at  Mg  =  4.8.  Here  we 
shall  present  analyses  of  the  recent  measurements  of  Holden  and  Hill.  Some 
preliminary  conparisons  were  given  in  Ref.  34,  but  the  diagnostic  techniques 
and  roughness  characterizations  have  recently  been  extended  and  refined,  per¬ 
mitting  more  definitive  analyses  at  this  time. 

One  interesting  series  of  tests  by  Holden®  were  performed  on  45° 
cones,  at  an  edge  Mach  number  of  1.8.  These  tests  illustrate  the  importance 
of  the  method  of  applying  grit  to  the  model  surface.  Holden's  most  recent 
studies  used  a  two-sided  tape  to  apply  the  particles,  whereas  earlier  methods 
used  Krylon  spray  adhesive.  Holden®  kindly  supplied  us  with  profilometer 
traces  of  the  surfaces  prepared  in  each  way  from  identical  "4  mil"  grit. 

Figure  9  shows  the  average  elements  derived  by  the  method  presented  in  Section 
3.  There  is  clearly  a  significant  difference,  with  the  two-sided  tape  yield¬ 
ing  greater  roughness  height  and  density.  Subsequent  discussions  with  Holden 
led  to  the  conclusion  that  various  effects  such  as  agglomeration  can  affect 
the  bonding  process.  It  is  not  sufficient  to  use  the  nominal  grit  or  sand 
grain  size  as  a  measure  of  the  roughness  size. 

In  Fig.  10  we  conpare  the  present  model  with  Holden's  measured  heat 
transfer  data  for  the  4-mil  tape  mounted  roughness.  A  significant  error  is 
seen  in  the  absolute  level  of  the  smooth  and  rough  wall  computations.  This 
error  seems  to  be  appreciable  for  cases  with  low  values  of  T^^/Tg,  at  modest 
edge  Mach  numbers.  At  hypersonic  speeds,  viscous  dissipation  maintains  high 
tenperatures .throughout  most  of  the  boundary  layer  and  no  such  overprediction 
is  found.  This  error,  then,  occurs  for  cases  with  extreme  density  or  tenpera- 
ture  variations  across  the  boundary  layer.  Despite  careful  excunination,  we 
have  not  been  able  to  remedy  this  discrepancy.  Numerical  accuracy  does  not 
seem  to  be  the  issue,  and  we  can  only  speculate  that  existing  turbulence  clo¬ 
sure  approximations  are  inadequate  in  situations  with  large  density  gradi- 
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If  one  were  to  scale  down  the  magnitude  of  the  calculations,  the  ex¬ 
tent  of  the  predicted  roughness  augmentation  would  agree  rather  well  with 
Holden's  data  for  spray-mounted  roughness.  As  also  indicated,  the  model  pre¬ 
dicts  less  increase  in  heating  for  the  Krylon-mounted  roughness. 

Two  very  significant  experiments  under  similar  hypersonic  conditions 
were  performed  by  Holden^^  and  Hill,^^  on  slender  cones  at  Mg  =  8-10.  Holden 
used  a  single  10  mil  (nominal)  grit,  emd  obtained  measurements  on  the  wind¬ 
ward  ray  at  angles  of  attac)c  from  0°  to  16°.  To  confirm  earlier  results  with 
thin  film  gages,  he  used  calorimeter  heat  transfer  gages,  in  addition  to  skin 
friction  gages.  Hill  used  three  different  roughnesses,  from  nominal  grit 
sizes  of  11,  37  and  65  mils.  Both  experimenters  provided  us  with  profilometer 
traces,  from  which  we  derived  the  roughness  parameters. 

Figure  11  shows  the  roughness  element  specifications  derived  from  the 
profilometer  traces.  Hill's  method  of  application  appears  to  result  in  ele¬ 
ments  that  are  more  vertically  aligned  and  more  closely  spaced  on  a  relative 
basis.  For  three  of  these  surfaces,  the  derived  roughness  height  is  close  to 
the  average  grit  diameter,  but  is  very  much  less  for  Hill's  "65  mil"  grit. 

This  finding  emphasizes  the  need  to  perform  careful  characterizations  of  actu¬ 
al  rough  surfaces. 

Figure  12  compares  tile  jresent  model  with  Holden's  measurements.  For 
the  cases  at  Mgle  of  attack,  we  used  an  equivalent  cone  approximation  to  de¬ 
scribe  angle  of  attack  effects,  which  might  be  less  accurate  at  higher  angles. 
Several  trends  are  evident  in  Fig.  12.  Roughness  causes  a  greater  increase  in 
skin  friction  than  in  heat  transfer.  At  a  =  0°,  Holden's  results  show  little 
roughness  effect  at  the  larger  downstream  distances,  whereas  a  modest  effect 
is  predicted.  Greater  increases,  and  better  agreement  between  theory  and 
data,  are  seen  with  increasing  angle  of  attack. 

The  corresponding  comparisons  with  Hill's  heat  transfer  data  are  shown 
in  Fig.  13.  The  predicted  roughness  heating  is  similar  for  all  three  of 
Hill's  roughness.  The  11  mil  case  is  somewhat  over-predicted,  while  the  other 
two  cases  are  underpredicted  at  greater  downstream  distances.  Both  theory  and 
data  indicate  slightly  higher  heating  values  at  37  mils  than  at  65  mils.  This 
behavior  is  compatible  with  the  derived  roughness  characteristics  shown  in 


Comparison  of  computed  skin  friction  versus  distance  with  Holden's  data'^'^  on  a  6’  cone 
at  M  =  9.4,  no  angle  of  attack. 
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Fig.  11,  which  showed  the  same  height  but  greater  spacing  for  the  65  mil  grit. 
One  puzzling  aspect  of  Hill's  data  is  the  tendency  of  the  rough  wall  heating 
rates  to  be  essentially  independent  of  distance.  If  anything,  we  would  expect 
the  rough  wall  measurements  to  decay  at  a  greater  rate  than  the  smooth  wall 
data,  since  k/B  decreases  with  increasing  distance. 

The  10  mil  data  of  Holden  at  0°  angle  of  attack  and  the  11  mil  data  of 
Hill  are  significant,  in  that  very  little  effect  of  roughness  is  evident  in 
either  case.  The  present  theory  and  most  existing  correlations  predict  at 
least  a  moderate  effect.  It  is  suggested  below  that  relatively  modest  rough¬ 
ness  Reynolds  numbers  are  responsible  for  this  behavior.  The  quantity  k"*"  is 
50-70  for  either  case,  but  is  much  larger  for  Holden's  cases  at  nonzero  angles 
of  attack  and  for  Hill's  cases  at  larger  roughness.  Such  k"*"  values  are  only 
slightly  below  the  fully-rcugli  requirement  of  k"*"  =  70  for  incon^jressible 
flows.  However,  careful  examination  of  our  conputer  results  indicates  that  the 
transition  values  should  increase  with  increasing  Mach  number,  to  be  discussed 
more  extensively  below. 


6.  ROUGHNESS  SCALING  LAW 


While  the  model  presented  above  yields  good  general  agreement  with  the 
available  data,  it  is  not  particularly  useful  for  engineering  purposes.  The 
conputational  cost  of  running  the  computer  code  is  modest.  However,  it  is  im¬ 
practical  to  expect  other  users  to  become  familiar  with  the  program,  which  in¬ 
volves  finite  difference  solution  of  many  simultaneous,  stiff  partial  differ¬ 
ential  equations.  What  is  needed  is  an  algebraic  recipe  that  can  be  readily 
understood  and  applied  to  practical  problems.  Our  computer  model  can  be  use¬ 
ful  in  developing  such  an  engineering  method  in  two  ways:  1)  numerical  solu¬ 
tions  can  be  examined  to  determine  the  dominant  physical  processes  and  2)  the 
code  can  be  exercised  to  generate  a  base  of  numerical  data  covering  the  range 
of  input  parameters  far  more  thoroughly  than  do  the  available  experiments. 

The  key  to  developing  scaling  laws  for  roughness  effects  lies  in  the 
confuted  mean  velocity  profiles.  For  the  vast  majority  of  cases  considered, 
the  mean  velocity  is  quite  uniform  over  much  of  the  range  y  <  k.  An  example 
is  shown  in  Fig.  14,  from  our  solution  for  Schlichting' s  closely  packed  spher¬ 
ical  segments.  Near  the  top  of  the  elements  the  velocity  profile  increases 

and  blends  into  the  log  region,  and  the  velocity  must  decrease  towards  zero  as 

y  ->■  0.  But  the  velocity  is  remarkably  constant  for  most  of  the  region  y  <  k. 
The  exceptional  cases,  in  which  a  range  of  uniform  velocity  is  not  evident, 
generally  involve  such  short  or  sparse  roughnesses  that  the  smoot.h  wall  pro¬ 
file  is  hardly  altered. 

It  must  be  admitted  that  this  velocity  behavior  was  not  anticipated  in 
our  model  development.  Physically,  for  y  <  k  the  turbulence  is  simply  diffus¬ 
ing  toward  the  base  of  the  wall  and  dissipating.  Perhaps  an  approximate  model 

could  be  developed,  for  example  by  a  three-layer  approach  (the  constant  veloc¬ 

ity  region,  the  log  region,  and  the  outer  region) .  There  is  little  experi¬ 
mental  evidence  to  confirm  or  deny  the  predicted  behavior,  since  it  is  gen¬ 
erally  not  feasible  to  measure  flow  properties  between  roughness  elements. 

Chen  and  Roberson show  a  profile  that  increases  by  only  10-20%  in  the  range 
0.2  <  y/k  <  1,  for  hemispheres  at  an  average  spacing  of  4.5  diameters.  This 
velocity  profile  behavior  is  apparently  well  known  in  the  study  of  turbulence 
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Fig.  14  Computed  mean  velocity  profile  for  a  typical  nse  (closely  packed  spherical  segments) . 


in  plant  canopies,  and  Raupach  and  Thoiti^^  quote  two  measured  wind  profiles  in 
plant  canopies  (a  pine  forest  and  a  maize  field)  that  are  quite  similar  to 
that  of  Fig.  14. 

The  wall  shear  is  given  by 


C 


f 


p  U  2 
e  e 


C  f(y) 
D 


(14) 


(the  bloc)cage  factor  f(y),  defined  in  Eq.  (10)  ,  enters  through  the  stream 
function) .  Now  let  us  approximate  U  as  a  constant  U  =  U|^.  For  conyressible 
cases,  p  will  be  related  to  U,  and  hence  is  also  constant  p  =  pj^.  The  block¬ 
age  factor  generally  varies  slowly  with  height  and  will  be  approximated  by  its 
value  at  k/2.  Since  the  integral  of  D(y)  is  simply  the  frontal  area  of  the 
roughness  element,  Eq.  (14)  reduces  to 
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which  is  the  basis  for  deriving  the  appropriate  scaling  laws. 

6.1  Inconpressible  Fully-Rough  Flow’ 

Equation  (15)  may  be  contrasted  to  Eq.  (2) 


(15) 


(2) 
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For  sand  grain  roughness/  the  velocity  shift  in  the  fully  rough  regime  is 

AU 

- - 5.6  log  V  -  3  . 


(3) 


The  smooth  wall  skin  friction  involves  the  Reynolds  number  based  on  some  meas¬ 
ure  of  the  width  of  turbulent  layer.  We  shall  select  the  momentum  thickness  6 
as  the  appropriate  thickness.*  To  a  very  good  approximation 


5.6  log  Reg  +  C  , 


(16) 


where  C  depends  on  flow  geometry 
(16)  into  Eq.  (2)  shows  that  the 
fully  rough  conditions 


and  pressure  gradient.  Putting  Eqs.  (3)  and 
skin  friction  is  solely  a  function  of  k/6  for 


5.6 


=  3  +  C  -  5.6  log  kg/e  (17) 

We  have  run  our  conputer  code  over  a  wide  range  of  parameters  such  ac 
roughness  shape  and  spacing.  We  first  confirmed  that  Cj  is  a  function  only  of 
k/6,  within  the  available  numerical  accuracy.  Figure  15  shows  an  exao^le  of 
the  conputed  behavior,  for  hemispherical  roughness  elements  at  three  relative 
spacings;  k"*"  varies  substantially.  The  spread  of  the  computed  points  for  each 
spacing  is  indicative  of  the  numerical  accuracy.  The  curves  represent  Eq. 
(17),  with  kg  chosen  to  match  the  computed  values  at  k/6  =  1. 


"The  boundary  layer  thickness  is  most  analogous  to  the  radius  of  a  pipe  or 
half-height  of  a  channel.  However,  boundary  layer  models  compute  6  inaccu¬ 
rately  (in  our  case)  or  not  at  all  (in  most  integral  methods).  The  displace¬ 
ment  thickness  6*  is  also  not  a  good  measure  for  compressible  situations,  for 
which  it  becomes  quite  sensitive  to  the  density  ratio  across  the  boundary 
layer. 
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Fig.  15 


Equations  (15)  and  (17)  represent  two  alternate  views  of  roughness  ef¬ 
fects.  Equation  (17)  explicitly  shows  the  dependence  on  kg/6,  but  contains  no 
information  on  roughness  character;  Eq.  (15)  is  useful  for  displaying  rough¬ 
ness  character.  If  Cf  is  to  depend  only  on  k/6  and  roughness  character/  then 
Ur  must  depend  on  the  same  quantities  and  must  also  depend  weakly  on  Re0  in 
such  a  fashion  as  to  cancel  the  weak  dependence  of  ^fgn,  on  Re0.  In  practice/ 
we  have  correlated  Ur  from  computed  values  of  C£  at  any  Re0  and  k"*"  in  the 
fully  rough  regime/  specifying  ^fgm  as  the  flat  plate  value  at  Re0  =  10^.  So 
long  as  the  same  value  of  ^fg®  is  retained,  the  resulting  correlations  for  Ur 
can  be  used  at  any  Re0  or  k"*"  in  the  fully  rough  regime. 

We  have  developed  a  numerical  data  base  for  roughness  character  ef¬ 
fects  from  boundary  layer  runs  with  a  variety  of  roughness  shapes  and  spac- 
ings.  The  shapes  considered  were  hemispheres,  spheres,  cylinders  (diameter 
=  height),  30°  (half-angle)  cones,  45°  cones,  and  truncated  30°  cones  (top 
diam  =  base  diam/2).  Spacing  was  generally  varied  from  k/l  =  0.1  to  as  close¬ 
ly  packed  as  numerically  feasible.  From  the  computed  values  of  Cf,  along  with 
Cq  =  0.6,  Pp/Pe  =  1,  ^fsm  “  1«81  ^  10“^  for  our  flat  plate  solution  at  Re0 
=  10^,  and  the  appropriate  values  of  f(k/2)  and  Xy.,  we  then  solved  Eq.  (15) 
for  UR/Ug.  An  example  of  the  derived  values  is  shown  in  Fig,  16  for  hemi¬ 
spheres.  The  straight  lines  are  drawn  as  an  aid  to  the  eye;  UR/Ug  cannot  be 
precisely  linear  in  log  k/0  and  still  have  consistency  between  Eqs.  (15)  f.nd 
(17.) 


The  other  shapes  investigated  yield  plots  quite  similar  to  that  of 
Fig.  16.  In  fact,  if  one  were  to  approximate  the  numerically  derived  values 
of  Ur  with  straight  lines  as  indicated  in  Fig.  16,  the  height  amd  slope  of  the 
lines  are  almost  solely  a  function  of  and  quite  insensitive  to  element 
shape.  Such  a  correlation  for  Ur  as  a  function  of  k/0  and  Xj^  could  then  be 
used  in  Eq.  (15)  as  a  useful  engineering  tool.  However,  it  is  more  precise 
and  more  compatible  with  existing  methods  to  correlate  Ur  vs.  X)^  at  a  fixed 
value  of  k/0  (such  as  unity)  and  then  use  Eqs.  (15)  and  (17)  to  derive  the  ef¬ 
fective  sand  grain  roughness.  In  Fig.  17  we  show  the" derived  values  of  Ur  at 
k/0  =  1  for  all  of  the  shapes  considered.  A  modest  shape  dependence  may  be 
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Fig.  16  Derived  values  of  roughness  velocity  versus  k/0,  for  hemispherical  roughness  at  various 


Fig.  17  Correlation  of  roughness  velocicy  at  k/8 


detected,  with  conical  elements  tending  to  fall  below  the  more  blunt  shapes. 
The  correlation  suggested  in  the  figure  best  fits  the  more  blunt  shapes,  which 
may  be  more  representative  of  actual  surface  roughnesses;  also,  our  theory 
tends  to  underpredict  the  roughness  effect  observed  for  cones  by  Schlicht- 


If  one  inserts  the  indicated  correlation 


~  =  0.247  +  0.234  log 
e 


(18) 


along  with  ^fsm  =  1.81  x  10“^,  C  =  5.24,  and  Cq  =  0.6  into  the  above  equa¬ 
tions,  the  following  equations  result  for  the  ratio  of  sand  grain  to  actual 
roughness  height. 


C  *  =  1.81  X  10-3  +  0.6(0.247  +  0.234  log  X)^)2  f(|)  -I”  {19a) 

f i  \2  I 

1/2 

5.6  log  kg/k  =  8.50  -  +  5.6  log^^^  (19b) 

Figure  18  compares  the  value  of  kg/k  predicted  by  Eq.  (19)  with  the 
actual  observed  values,  for  a  number  of  experiments  involving  roughness  char¬ 
acter.  Two  sets  of  data  are  seriously  underpredicted  by  the  present  correla¬ 
tion  (as  well  as  by  the  computer  model  -  cf.  Figs.  6  and  8) s  Schlichting' s 
short  angles^**  (baffles)  and  the  stones  of  Mirajgoaker  and  Charlu. We  have 
no  explanation  for  these  discrepancies.  The  tall  rods  of  Seginer  et  al.22  and 
Thora21  are  rather  substantially  underpredicted. 


For  comparison.  Fig.  19  shows  a  similar  comparison  for  the  correlation 
of  Simpson, 26  a  refinement  of  Dvorak’s  correlation. 2^  The  general  degree  of 
correlation  is  comparable  to  that  of  the  present  result.  The  short  angles  of 
Schlichting  and  stones  of  Mirajgoaher  and  Charlu  are  again  underpredicted;  the 
rods  of  Seginer  are  handled  decently,  but  those  of  Thom  could  not  even  be 
plotted  on  the  figure.  It  is  roost  important  to  note  that  the  Simpson  (and 
Dvorak)  correlation  is  not  accurate  for  closely  packed  roughness  elements. 
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Fig.  19  Comparison  of  predicted  versus  actual  roughness  for 
Simpson’s  correlation. 


The  points  for  Schlichting* s  closely  packed  spheres,  spheres  with  l/V  *  1.5, 
and  closely  packed  spherical  segments  all  fall  outside  the  indicated  ±  20% 
error  band  in  Fig.  19.  Most  of  the  available  experimental  data  on  roughness 
character  have  been  obtained  with  relatively  sparse  roughness  patterns.  The 
Simpson/Dvorak  correlations  are  fine  tuned  to  the  data  base,  but  do  not  have 
the  benefit  of  sufficient  distributed  roughness  data  at  greater  roughness  den¬ 
sities.  Use  of  the  straight  line  segment  with  positive  slope  in  Fig.  1  leads 
to  the  poor  results  cited  and  Simpson  questioned  the  validity  of  that  portion 
of  the  correlation  for  3-D  roughness.  Hence,  for  relatively  closely-packed 
roughness  (say  i/k  <  1.5),  which  corresponds  to  most  surfaces  of  practical  in¬ 
terest,  the  correlation  derived  here  should  be  the  more  reliable. 

6.2  Compressible  Flows-Fully  Rough 

Compressibility  effects  generally  reduce  friction  on  both  smooth  eoid 
rough  walls.  Smooth  wall  friction  coefficients  are  generally  computed  from 
"transformation  functions,"  Fg  and  F^.  As  reviewed  in  Ref.  29,  an  equivalent 
incompressible  coefficient  is  computed  at  the  effective  Reynolds  number  F0  • 
Re9,  and  is  then  divided  by  F^  to  obtain  the  compressible  coefficient: 

^fsm  =  f-  C  (Fe  •  Ree)  (20) 

'c  *i'®m 

As  noted  above,  Dvorak^‘*»25  used  Eq.  (2)  to  obtain  a  rough  wall  incom¬ 
pressible  friction  coefficient,  evaluating  ^fgm  ^6  *  *^^6  AOi/Ux  at  the 
actual  value  of  k"*";  he  then  divides  by  Fj,  according  to  Eq.  (20).  Note  that  Fg 
plays  a  very  small  role  -  F0  is  usually  near  unity  and  ^fgn,  is  insensitive  to 
Reynolds  number.  If  Fg  plays  no  role,  then  this  procedure  guarantees  that  the 
percentage  increase  of  friction  is  independent  of  compressibility,  at  fixed 
k/0.  Our  model,  however,  predicts  a  rather  different  result  -  that  the  rela¬ 
tive  effect  of  roughness  decreases  with  increasing  Mach  number. 

In  the  present  model,  one  obvious  effect  of  compressibility  is  through 
the  density.  At  high  Mach  numbers,  viscous  dissipation  causes  high  tempera¬ 
tures  within  the  Ixjundary  layer.  Because  roughness  reduces  velocities  within 


the  boundary  layer,  even  higher  temperatures  can  be  expected.  To  estimate  the 
roughness  density,  we  evaluated  the  relation  between  total  enthalpy  and  veloc¬ 
ity  in  the  output  from  a  large  number  of  computer  runs  at  edge  Mach  numbers 
between  0  and  10.  A  linear  relationship  was  found  to  be  quite  accurate 
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(21) 


as  it  is  with  nearly  any  turbulent  situation.  Then,  with  perfect  gas  rela¬ 
tions  it  follows  easily  that 


e 


Figure  20  compares  this  equation  with  the  values  of  pj^  (at  y  =  k/2)  and  Uj^ 
from  several  runs  of  the  computer  model.  Except  at  Mg  =  10,  T^/Tg  =  0.74, 
where  Eq.  (22)  is  perhaps  20%  high,  the  equation  provides  a  simple  and  accu¬ 
rate  representation  of  the  density. 

The  roughness  density  term  can  obviously  significantly  reduce  the  mag¬ 
nitude  of  the  friction  predicted  by  Eq.  (15).  This  reduction  is  similar  to, 
but  not  precisely  the  same  as,  that  obtained  by  dividing  an  incompressible 
value  by  Fg.  To  a  good  approximation  (as  in  the  Sommer-Short  method^^),  Fg  is 
given  by  the  reference  temperature 

T  T 

rftf  w  2 

F  »  =  0.55  +  0.45  +  0.035  M  .  (23) 

c  T  T  e 

e  e 


-49- 


4, 


The  value  of  Tj-gf  tends  to  be  somewhat  below  the  values  of  Tj^  given  by  Eq. 

(22)  or  Fig.  20  -  for  T^/Tg  =  0.74  they  agree  at  UR/Ug  =  0.23;  at  T^^/Tg  = 
they  agree  at  Uj^/Ug  =  0.12.  However,  the  computer  model  indicates  an  addi¬ 
tional  effect,  in  that  U|^  is  affected  by  compressibility. 

A  careful  examination  of  the  values  of  UR/Ug  derived  from  our  computer 
solutions,  over  a  wide  range  of  Mach  numbers  and/or  wall  temperature  ratios, 
indicates  that  the  roughness  velocity  scales  with  PrIjc/Pb/  indicated  in 
Fig.  21.  Note  that  most  of  the  scatter  about  the  incompressible  correlation 
results  from  cases  where  the  wall  temperature  ratio  was  varied;  our  computer 
code  has  inherent  inaccuracies  for  cases  with  T^,/Tg  very  small  or  large,  even 
for  smooth  walls.  Otherwise,  Fig.  21  shows  a  rather  solid  correlation  for  the 
effect  of  compressibility  on  Uj^.  Note,  however,  that  this  correlation 

—  =  0.247  +  0.234  logl  —  j  (24) 

e  '  *^e 


is  a  transcendental  equation  for  Ur,  since  Pr  depends  on  Ur  through  Eq.  (22). 
But,  Pr  is  typically  insensitive  to  Ur,  and  one  can  iterate  to  a  solution  very 
quickly. 

Computation  ol  the  compressible  skin  frictioij  requires  first  the  com¬ 
pressible  coefficient  at  k/0  =  1,  analogous  to  the  incompressible  value  given 
by  Eq.  ( 19a) : 


Fig.  21  Correlation  of  roughness  velocity  at  k/0  =  1  for  compressible  cases. 
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Again,  we  have  a  transcendental  equation,  but  one  that  is  easily  solved  itera¬ 
tively. 

If  one  prefers  to  express  the  emswer  in  terms  of  an  equivalent  sand- 
grain  roughness  (from  incompressible  data  or  Eqs.  (19)},  then  there  must  be  an 
upward  shift  in  the  velocity 


(^)  -it) 


Here,  ^fsm  is  the  actual  compressible  smooth  wall  value  and  AU^/U^  is  calcu¬ 
lated  from  the  low  speed  relation  (Eq.  3)  using  the  actual  wall  conditions  for 
)«■'■.  The  "compressibility  shift"  depends  on  roughness  character  and  compress¬ 
ibility  conditions  (M^,  T^/Tg): 
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_  _  /f^  V  \ 


-  8.5  /f~  +  5.6/F  ~ 

c  c  ^VWp  V 

\’  w  w 


where  Cf  comes  from  Eq.  (25). 
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This  shift  can  be  substantial  -  for  example,  for  hemispheres  at  k/f 
=  0.4,  Mg  =  10,  T„/Tg  =  0.74,  the  value  is  14.3. 

6.3  Rough/Smooth  Transition 

With  Nikuradse's  sand  grain  roughness,  the  transition  between  rough 
and  smooth  wall  behavior  occurs  in  the  range  5  <  kg'**  <70.  It  is  questionable 
whether  the  sand  grain  behavior  could  be  applied  to  cases  with  varying  rough¬ 
ness  character  or  Mach  number.  Dvorak^^  presents  an  interpolation  method  for 
the  transitional  regime,  but  there  is  very  little  supporting  data  for  3-D 
roughness.  We  have  investigated  this  issue  with  our  computer  model,  recogniz¬ 
ing  the  likely  limitations  of  the  theory.  The  form  drag  assumption  is  basic¬ 
ally  a  high  Reynolds  number  concept.  Viscous  drag  on  the  elements  is  simply 
lumped  into  the  underlying  smooth  wall  friction,  and  there  is.no  reason  to  be¬ 
lieve  that  our  computer  model  accurately  describes  drag  on  roughness  elements 
at  lower  Reynolds  numbers.  Interference  between  neighboring  elements  also  be¬ 
comes  more  important  with  decreasing  Reynolds  numbers.  However,  the  results 
are  certainly  interesting. 

Figure  22  shows  the  values  of  Nikuradse's  quantity  B  (-3  is  replaced 
by  5.5-B  in  Eq.  (3))  derived  from  our  computer  calculations  for  hemispherical 
roughness  at  three  spacings.  The  derived  values  depart  from  the  fully  rough 
behavior  (B  =  8.5)  at  approximately  the  same  value  of  kg'*"  for  all  three  spac¬ 
ings.  The  more  dense  patterns  show  a  transitional  behavior  similar  to  that 
observed  by  Nikuradse,  at  least  within  the  accuracy  of  our  calculations.  How¬ 
ever,  the  solutions  for  wide  spacing  show  a  completely  different  behavior. 

The  transitional  behavior  is  even  more  complex  for  supersonic  flows. 
Our  solutions  consistently  indicate  that  minimum  k"*"  value  for  fully  rough  be¬ 
havior  increases  with  increasing  Mach  number.  In  Fig.  23  we  show  curve  fits 
through  the  derived  values  of  B  with  increasing  edge  Mach  number,  for  hemi¬ 
spheres  at  k/l  =>  0.4  (the  scatter  of  the  computer  data  is  substantial,  but  the 
trends  are  clear).  Two  effects  are  apparently  involved  here.  First,  the 
smooth  wall  solution  tends  to  shift  to  increasing  valQes  of  k'*'  with  increasing 
Mach  number,  as  the  compressibility  term  increases.  Second,  as  dis¬ 

cussed  edx3ve  with  regard  to  fully  rough  behavior,  the  effective  temperature  in 


Fig.  22  Computed  rough-smooth  transition  behavior  in  low  speed  flow. 


the  vicinity  of  the  roughness  elements  should  reflect  the  substantial  viscous 
dissipation  that  occurs  at  high  edge  Mach  numbers.  Even  for  relatively  low 
values  of  k"*",  according  to  our  solutions,  the  fluid  properties  at  y  =  k/2  or 
y  =  k  will  differ  significantly  from  wall  conditions,  upon  which  k'*’  is  based. 
For  example,  for  the  Mach  10  case  shown  in  Fig,  23  at  kg"*"  =  300,  where  the 
solution  shows  definite  departure  from  the  fully  rough  regime,  the  temperature 
at  y  =  k/2  is  about  3.6  Tg  or  4.8  With  properties  based  on  this  temper¬ 

ature,  rather  than  T^,  the  resulting  value  of  kg"*"  would  be  reduced  by  a  factor 
of  6,  bringing  it  into  reasonable  alignment  with  the  low  speed  curve. 

Figure  24  shows  the  computed  velocity  shifw  for  conditions  correspond¬ 
ing  to  Holden's  experiments  at  zero  angle  of  attack,  plotted  against  k"*"  (k 
=  10  mil,  kg  =  22  mil  from  Eq.  (19b)).  As  indicated,  the  test  conditions  cor¬ 
respond  to  k"*"  -  50-70,  where  the  solution  definitely  departs  fiom  the  fully 
rough  solution.  The  computer  model  yields  a  value  about  35%  below  the  ful¬ 
ly  rough  value  for  the  same  k/6.  The  fact  that  Holden’s  measurements  are  vir¬ 
tually  indistinguishable  from  smooth  wall  values  suggests  that  the  true  depar- 
txire  from  fully  rough  behavior  is  even  greater  than  predicted  at  k'^  =  50-70. 
The  conditions  for  Hill's  11  mil  roughness  are  very  much  the  same.  Thus,  we 
conclude  that  a  combination  of  compressibility  and  transitional  (smooth/rough) 
effects  combine  to  cavise  the  minimal  augmentation  observed  in  both  cases. 
However,  it  would  be  foolhardy  to  attempt  to  derive  a  detailed  description  of 
the  rough/smooth  transition  zone  from  the  present  model,  and  a  much  better 
data  base  is  clearly  needed. 

6.4  Rough  Wall  Prandtl  Number 

It  is  commonly  observed  that  roughness  causes  a  smaller  augmentation 
of  heat  transfer  than  of  skin  friction.  In  terms  of  the  present  model,  this 
results  from  the  fact  that  there  is  no  thermal  analogy  to  form  drag.  We  pre¬ 
viously^^  noted  that  all  components  of  the  fluctuating  velocity  are  in¬ 
creased  proportionally  by  roughness,  while  the  fluctuating  temperature  is  es¬ 
sentially  unchanged  from  the  smooth  wall  value.  This  reasoning  suggests  that 
the  augmentation  in  the  heat  flux  (-v’h')  is  the  square  root  of  that  of  the 
friction  (-u'v*) 
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Alternatively,  the  empirical  result  of  Dahm  et  al.^®  is 
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In  Fig.  25  we  show  the  computed  Stanton  numbers  and  skin  friction  co¬ 
efficients  for  various  cases  spanning  variations  in  roughness  height,  rough¬ 
ness  shape  and  Mach  number.  The  scatter  of  computer  points  may  be  partially 
indicative  of  the  inherent  numerical  accuracy  of  our  computer  model,  although 
the  points  showing  the  largest  departure  from  the  mean  generally  correspond  to 
upstream  locations,  where  the  solution  still  reflects  initial  conditions. 


The  Dahm  result,  Eq.  (30),  is  seen  to  provide  an  excellent  fit  to  the 
computer  solutions,  particularly  at  small  augmentation  ratios.  The  square 
root  law  of  Eq.  (29)  has  an  appropriate  functional  dependence,  but  consistent¬ 
ly  underpredicts  the  computed  heating  augmentation.  A  better  curve  fit  to  the 
various  computed  points  shown  in  Fig.  25  is  a  combination  of  the  functional 
dependence  of  Eqs.  (29)  and  (30): 
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Fig.  25  Computed  l.^at  transfer  augmentation  versus  skin  friction  augmentation. 
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X  [Roughness  base  area/unit  area]“^ 

Xj^  [Roughness  frontal  area/unit  area]~^ 

y  Dynamic  viscosity 

V  Kinematic  viscosity/  p/p 
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